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Abstract 

As of version 8.150 of Pythia, the isotropic decay model of r-leptons has been replaced with sophisticated r-lepton 
decay machinery. The decays and spin correlations for r-leptons in Pythia 8 are described, including the spin corre- 
lation algorithm, the available r-lepton production processes, the r-lepton decay models, the user interface, and the 
implementation. 
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1. Introduction 

The role of r-leptons in Higgs boson measure- 
ments [1 1 and beyond the Standard Model searches [2 | 
is becoming increasingly important, due to the en- 
hanced coupling of the r-lepton in many of these 
physics models. Consequently, it is necessary for cur- 
rent Monte Carlo generators to ensure accurate model- 
ing of r-lepton decays. 

Prior to version 8.150 of Pythia 0, r-lepton decays 
in Pythia were performed using a leptonic or generic 
hadronic current matrix element without including spin 
correlations, and more sophisticated r-lepton modeling 
was only possible through external packages such as 
Tauola |4|. Now, in Pythia version 8.150 and above, 
fully modeled hadronic currents with spin correlations 
are available, based on prior r-lepton modeling work 
in Tauola and Herwig++ 0. Currently, all known 
r-lepton decays with a branching fraction greater than 
0.04% are modeled. 

In this review, the spin correlation algorithm used in 
Pythia for r-leptons is described, as well as the avail- 
able r-lepton production mechanisms, summarized in 
Table [T] and r-lepton decays, given in Table [2] An in- 
troduction to the user interface, as well as the technical 
implementation is given. 

2. Correlations 

The spin correlation algorithm used in Pythia for r- 
lepton decays is based on the algorithm proposed by 
Collins [6 1 and Knowles [7|, and expanded by Richard- 
son J8). The algorithm separates spin correlations from 



the hard process, parton shower, and hadronization 
phases of the Monte Carlo generator, while maintaining 
full correlations, and can be broken into the following 
steps. 

1. The 2 — » n hard process is generated according to 
its matrix element At 

2. One of the outgoing particles is selected and its he- 
licity density matrix is calculated, 

*rte« m 

and the trace is normalized. Here, p 1 ' 2 are the helic- 
ity density matrices of the incoming particles with 
helicity k\ M is the matrix element with outgo- 
ing particle helicities Ak, and D k are the decay ma- 
trices of the outgoing particles, initialized to the 
identity. If the particle is from a decay, the helicity 
density matrix is the same as above but without p 2 
and K2- 

3. The selected particle is decayed using the weight, 

k=\,n 

where p is the helicity density matrix of the decay- 
ing particle with helicity Ao and the decay matrix 
element M. 

4. Steps [2] through [3] are performed until a decay is 
reached with no unstable particles. 



Preprint submitted to Nuclear Physics B Proceedings 



November 29, 2012 



c 

o 



c 

3 



C 

X3 



o 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 






■ ■i 


■ v s„ 

- "m 


■Z V- ! 

J ! 

// r 

.... W J* . 




i 




r ^ 

- -l.-1._j 


b ./" 




i / 




p- i i i I i i i I 


>J 



0.2 0.4 0.6 0.8 1 

2E K -I^ 

Figure 1 : Distribution of the fractional energy of the pion from a r~ — » 
v T 7r~ decay in the rest frame of the producing boson. The r-leptons 
are produced from Z (dashed black), H (dotted red), W (dash-dotted 
green), and H ± (dash-dot-dotted blue) bosons. 



Because the Higgs is spin zero, the production in- 
formation is not necessary for calculating r-lepton spin 
correlations. However, the type of Higgs must be 
known. Currently, spin correlations of r-leptons pro- 
duced from the CP-even (H, h°, H°), CP-odd (A ), or 
charged (if*) Higgses are calculated. 

The r-lepton can also be produced from a variety of B 
and D meson decays. For these decays, the constituent 
quarks of the meson are approximated, and the spin cor- 
relations are calculated using an //' — > W — > //' ma- 
trix element. For any r-lepton from an unknown pro- 
cess, the r-lepton is assumed to be unpolarized. 

The slope of the energy distribution for pions from 
r~ — > v T 7r~ decays in the rest frame of the parent 
bosons is proportional to the average r-lepton polariza- 
tion. In Figure [T] this distribution, from Pythia simula- 
tion, is given for r-leptons produced from electroweak 
and Higgs processes. Note the the opposite polarization 
of r-leptons produced from a W and H ± , the unpolar- 
ized r-leptons produced from neutral Higgses, and the 
slight polarization of r-leptons produced from Z bosons 
at a pp collider. 



5. The decay matrix of the last decayed particle is cal- 
culated, 

k=l,n 

and the trace is normalized. 

6. An undecayed particle from the decay above is ran- 
domly selected and steps [2] through [5] are repeated. 

7. Step [6] is repeated until all unstable particles are 
decayed. 



Type 



Processes 



electroweak 



// -» y -> //. // -> z -> //, 

// -> y'lZ -> //, //' -> W -> //', 
Z -» //, W -> ff 



Higgs 



H -> ff, h° -> ff, H° 
A -. //, H+- -» ff 



ff 



other 



B/D -> ff + X 



Table 1 : Production mechanisms in Pythia for which full spin corre- 
lations of r-lepton decays are automatically calculated. 



3. Production 

Spin correlated r-lepton decays can be generated in 
Pythia from a variety of hard processes, summarized in 
Table [T] The spin of the r-lepton can also be provided 
from an externally generated hard process, or explicitly 
set for all r-leptons or r-leptons from a specified parent. 

For electroweak processes, spin correlations are fully 
handled if the incoming fermions producing the y, Z, or 
W are known. If the W is not produced from fermions 
or the fermions are unavailable, the r-lepton is given a 
polarization of P = — 1, e.g. H — > W — > rv T . For y 
or Z hard processes not produced from fermions or with 
unavailable fermion information, the y or Z is assumed 
to be unpolarized, e.g. H — > ZZ — > 4r. 



4. Decays 

The matrix element for the decay of the r-lepton can 
be written as, 

M = /|- V (4) 

where gw is the SU(2) coupling, mw the W mass, L M 
the leptonic current of the r-lepton, and a leptonic 
or hadronic current dependent upon the decay. The r- 
lepton current is S Vr 7 /J (l - y 5 )u T where u T is the r- 
lepton spinor, dependent upon momentum and helicity, 
and u Vr is the r-lepton neutrino spinor. A full list of 
the available currents modeled in Pythia is given in 
Tabled 
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For the two-body decays of the r-lepton, is fq, 
where / is a constant and q is the momentum of the 
hadron. For the three-body leptonic decays, r~ — > 
v T e~v e and r~ — > v T pTv^ is of the same form as L„. 
Two hadronic three-body decay models are available, a 
decay via a vector resonance and a decay via a vector 
and scalar resonance. 

Four-body r-lepton decays in Pythia are imple- 
mented in four different models. The primary four-body 
decays, r~ — > v T 7r°7r°7!-~ and r~ — > v T n~n~n + , are mod- 
eled using the CLEO fit. These decays can also be per- 
formed using a generic hadronic four-body model from 
Decker, et al. The four-body decays with kaons in the 
final state are calculated using a model from Finkemeier 
and Mirkes. 

The five-body decays of the r-lepton to pions are pro- 
duced with the Novosibirsk model, a phenomenological 
fit of four pion production from electron-positron anni- 
hilation. The six-body decays of the r-lepton are han- 
dled through a model proposed by Kiihn and Was. 

A comparison of the hadronic invariant mass distri- 
butions, generated with Pythia, between the dominant 
three, four, five, and six-body decays of the r-lepton is 
given in Figure [2] The r~ — > v T n°7t°7T~ distribution was 
generated using the default CLEO model. The expected 
upwards shift for higher multiplicity decays can be seen, 
as well the distinct difference in shapes due to the dif- 
fering propagators of the models. 



5. Interface 

Documentation for r-lepton decays in Pythia can 
be found under the Tau Decays subsection of Particle 
Decays in the Pythia HTML manual The r-lepton 
decay mechanism is configured by setting the mode 
sophisticatedTau of ParticleDecays. This mode 
default is 1, where the decays are performed using the 
sophisticated r-lepton decay machinery. When sophis- 
ticated decays are not necessary, or a comparison with 
the isotropic Pythia model is needed, the mode can be 
set to 0. 

If the mode is set to 1, and a LHEF file is read, the 
spin of the r-lepton will be set to the SPINUP informa- 
tion. However, if no SPINUP information is provided, 
and the r-lepton is from a known production process, 
Pythia will automatically calculate the polarization. If 
the r-lepton is from an unknown process the sophisti- 
cated decay will be performed assuming the r-lepton is 
unpolarized. 

It is possible to force the polarization of the r-leptons 
and override either the Pythia calculated polarization or 
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Figure 2: Distribution of the invariant mass of the hadronic current 
for the dominant three-body n°7T~ (dashed black), four-body n°n°7T~ 
(dotted red), five-body n°n~n~7i + (dash-dotted green), and six-body 
n~TCiC-n*n + (dash-dot-dotted blue) decays of the r-lepton. 



the polarization from a LHEF file. If the mode is set 
as 2, all r-leptons produced from a parent with an iden- 
tification code specified by tauMother will be given a 
polarization set by the parameter tauPolarization. If 
mode 3 is selected, all r-leptons will be produced with 
the polarization set by tauPolarization. 

The models used to decay the r-leptons for each 
channel can be changed by switching the matrix el- 
ement mode of the channel. The syntax takes the 
form 15 : channel :meMode = mode where 15 is the r- 
lepton particle identification code, channel specifies 
the decay channel number as listed under Particle Data, 
and mode is the new matrix element mode. For example, 
the default CLEO model used for the r~ — > v T 7r 7r°7r~ 
decay can switched to the Decker, et al. model using 
15:9:meMode = 1543. However, caution is advised 
in switching models, as it is possible to set a channel 
to use an incompatible matrix element, causing unpre- 
dicted behavior. 



6. Implementation 

The implementation of the r-lepton decay machinery 
in Pythia is intended to be complete and fast, yet eas- 
ily extensible for new decay models or production pro- 
cesses. A helicity framework is available in the Pythia 
source files HelicityBasics where methods for cal- 



3 



Multiplicity 


Model 




M 


Decay Products 


2-body 


single hadron 




1521 


n-,K- 




leptonic 




1531 




3 -body 


Kiihn and Santamaria 


® 


1532 


7T°7T-, K°K-, T]K- 




Pinkemeier and Adirkes 




1 JJJ 


71 A. , 71 A. 




CLEO 


El 


1541 


7T°7T°7T~, 7T~7r~7T + 


4-body 


Finkemeier and Mirkes 
Decker, et al. 


E3 
ED 


1542 
1543 


K-n-K\ K°n-K°, A^tTA^, A^ttA^, K°7t-K l , 
K-7r°K Q , n°n°K-, K-nn + \ n-R°7x 
n°n°n + , n-n-n + , K n K + , A°7r A , K-n°K°, 
AT, K-n-7T + , 7rAV, n-^T} 




Jadach, et al. 


03 


1544 


y7i°7i~ 


5 -body 


Novosibirsk 


E2 


1551 


7T 7T-n-n + ,n n°n n- 


6-body 


Kiihn and Was 


ED 


1561 


7l°7T°7T-7l-7T + , AV/^, n H n V TT* 



Table 2: Summary of available r-lepton decay models in Pythia 8 sorted by multiplicity. For each model the internal Pythia matrix element mode 
identifier (M) is given, as well as the decays available through the model and the reference of the model. The implicit v T is omitted for brevity. 



culating the helicity dependent wavefunctions for parti- 
cles, both fermions and bosons, is provided, as well as 
the corresponding y-matrices. 

Using this framework, the helicity matrix element for 
the decay of a r-lepton into a pion, 

M = iiy^il - y 5 )u T f n cf„ (5) 

can be written in pseudo-code as, 

M = ^v T .waveBar(/l Vr ) * GammaMatrix(/z) 

^ *(1 - GammaMatrix(5))*r.wave(/} T )*/ ff (6) 
*GammaMatrix(4)(//,^) * Wave4(g jr )( J u) 

where y 4 provides the H metric. 

Both the production matrix elements and r-lepton de- 
cay matrix elements are written using this framework 
and are provided in the HelicityMatrixElements 
files of the Pythia source. For the r-lepton decay matrix 
elements, it is only necessary to implement the cur- 
rent. The spin correlation and decays of the r-leptons 
are implemented in the TauDecays files of the Pythia 
source. 

7. Conclusion 

Fully modeled r-lepton decays with spin correla- 
tions are now available in Pythia, using fully extensi- 
ble code. Further development of r-lepton production 
mechanisms and decay models is underway, and feed- 
back on requested features is welcome. 
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